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We analyze the entanglement degradation provoked by the Hawking effect in a bipartite system 
Alice-Rob when Rob is in the proximities of a Schwarzschild black hole while Alice is free-falling into 
it. We will obtain the limit in which the tools imported from the Unruh entanglement degradation 
phenomenon can be used properly, keeping control on the approximation. As a result, we will be 
able to determine the degree of entanglement as a function of the distance of Rob to the event 
horizon, the mass of the black hole, and the frequency of Rob's entangled modes. By means of this 
analysis we will show that all the interesting phenomena occur in the vicinity of the event horizon 
and that the presence of event horizons do not effectively degrade the entanglement when Rob is 
far off the black hole. The universality of the phenomenon is presented: There are not fundamental 
differences for different masses when working in the natural unit system adapted to each black hole. 
We also discuss some aspects of the localization of Alice and Rob states. All this study is done 
without using the single mode approximation. 

PACS numbers: 03.67.Mn, 03.65.-w, 03.65.Yz, 04.62. +v 



I. INTRODUCTION 

Rclativistic quantum information is born from the 
combination of very different and fruitful branches of 
physics. Namely, general relativity, quantum field the- 
ory and quantum information theory. Its aim is to study 
the behavior of quantum correlations in rclativistic set- 
tings. In its scope, among other topics, is the study of the 
behavior of quantum correlations in non-inertial settings, 
which has produced abundant literature [H4l9| . This dis- 
cipline provides novel tools for the analysis of the Unruh 
effect and the Hawking effect, allowing us to study the 
behavior of the correlations shared between non-inertial 
observers. 

In previous works it was analyzed the entanglement 
degradation phenomenon produced when one of the part- 
ners of an entangled bipartite system undergoes a con- 
stant acceleration; this phenomenon, sometimes called 
Unruh decoherence, is strongly related to the Unruh ef- 
fect. Its study revealed that there are very strong differ- 
ences between fermionic and bosonic field entanglement 
d, [H, [13, Eil ■ The reason for these differences was 
traced back to fermionic/bosonic statistics and not to the 
difference between bosonic and fermionic mode popula- 
tion as previously thought (20l - [22j . In these earlier stud- 
ies some conclusions were drawn about the infinite accel- 
eration limit, in which the situation is similar to being 
arbitrarily close to an event horizon of a Schwarzschild 
black hole. 

However there arc many subtleties and differences be- 
tween Rindler and Schwarzschild space-times. For exam- 
ple Schwarzschild space-time presents a real curvature 
singularity while Rindler metric is nothing but the usual 
Minkowski metric represented in different coordinates 
and, therefore, has no singularities. The Rindler horizon 
is also of very different nature from the Schwarzschild's 



event horizon. Namely, the Rindler horizon is an acceler- 
ation horizon experienced only by accelerated observers 
(at rest in Rindler coordinates). On the other hand, a 
Schwarzschild horizon is an event horizon, which affects 
the global causal structure of the whole space-time, in- 
dependently of the observer. Also, for the Rindler space- 
time there are two well defined timelike Killing vectors 
with respect to which modes can be classified according 
to the criterion of being of positive or negative frequency. 
Contrarily, Schwarzschild space-time has only one time- 
like Killing vector (outside the horizon). 

Therefore, to analyze the entanglement degradation 
produced due to the Hawking effect near a Schwarzschild 
black hole we must be careful, above all if we want to 
do a deeper study than simply taking the limit in which 
the Rindler acceleration parameter becomes infinite. In 
this paper we will show how we can use the tools coming 
from the study of the Unruh degradation in uniformly ac- 
celerated scenarios without restricting only to the exact 
infinite acceleration limit and controlling to what extent 
such tools are valid. 

Consequently, we will be able to compute the entan- 
glement degradation introduced by the Hawking effect as 
a precise function of three physical parameters, the dis- 
tance of Rob to the event horizon, the mass of the black 
hole, and the frequency of the mode that Rob has en- 
tangled with Alice's field state. As a result of this study 
we will obtain not only the explicit form of the quan- 
tum correlations as a function of the physical parameters 
mentioned above, but also a quantitative control on what 
distances from the horizon can be still analyzed using the 
mathematical toolbox coming from the Rindler results. 

Contrarily to all the previous works in which it was 
carried out the single mode approximation described in 
0, [l2[ and of common use in the published literature, we 
will argue that we do not need to use such approximation 
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to study the fundamental effects on the entanglement 
derived from the Hawking effect. 

Our setting consists in two observers (Alice and Rob), 
one of them free-falling into a Schwarzschild black hole 
close to the horizon (Alice) and the other one standing at 
a small distance from the event horizon (Rob) . Alice and 
Rob are the observers of a bipartite quantum state which 
is maximally entangled for the observer in free fall. The 
Hawking effect will introduce degradation in the state as 
seen by Rob, impeding all the quantum information tasks 
between both observers. 

In this context we will analyze not only the classical 
and quantum correlations between Alice and Rob, but 
also what is the behavior of the correlations that both 
observers would acquire with the mode fields on the part 
of the space-time that is classically unaccessible due to 
the presence of the event horizon. 

By means of this study we will show that all the in- 
teresting entanglement behavior occurs in the vicinity of 
the event horizon. What is more, we will argue that as 
the entangled partners go away from the horizon the ef- 
fects on entanglement become unnoticcably small and, as 
a consequence, quantum information tasks in universes 
that contain event horizons are not jeopardized. 

We will also show that the phenomenon of the Hawk- 
ing degradation is universal for every Schwarzschild black 
hole, which is to say, it is ruled by the presence of the 
event horizon and is not fundamentally influenced by the 
specific value of the black hole parameters when the anal- 
ysis is performed using natural units to the black hole. 
Furthermore, we will discuss the validity of the results 
obtained when instead of the usual plane wave basis we 
work in a base of wave packets, for which the states of 
Alice and Rob can be spatially localized. 

This paper is structured as follows. In section [H] we 
show how we work without using the single mode approx- 
imation, presenting previous results about the Unruh en- 
tanglement degradation for scalar and Dirac fields. In 
section Hill we study the entanglement in a Schwarzschild 
space-time using the tools built for the Rindler case, de- 
tailing to what extent this approximation holds. In sec- 
tion IIVI we will present the result for the correlations 
between the different bipartitions in the Schwarzschild 
space-time scenario. In section [V] we show that the re- 
sults obtained in the preceding sections are also valid 
when we consider complete sets of localized modes in- 
stead of plane waves bases. Finally, we present our con- 
clusions in section |VT1 



II. RINDLER SPACE-TIME 

Along this work we are going to consider bipartite 
scalar and Dirac field states. We will name Alice the 
observer of the first part of the system and Rob the ob- 
server of the second part. In this fashion, the quantum 
state for the whole system is defined by the tensor prod- 
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Figure 1: Flat space-time conformal diagram showing Alice, 
Rob and AntiRob trajectories. i° denotes the spatial infini- 
ties, i~ , i + are respectively the timelike past and future in- 
finities, J ! ~ and J'^ are the null past and future infinities 
respectively, and are the Rindler horizons. 

uct 

|^A,<M = |^a) ®|^r). (1) 
Alice's Rob's 

Now, while Alice is in an inertial frame, we will consider 
that Rob is observing the system from an accelerated 
frame. 

A uniformly accelerated observer viewpoint is de- 
scribed by means of the well-known Rindler coordinates 
[23l [H| ■ In order to map field states in Minkowski space- 
time to Rindler coordinates, two different sets of coordi- 
nates are necessary. These sets of coordinates define two 
causally disconnected regions in Rindler space-time. If 
we consider that the uniform acceleration a lies on the z 
axis, the new Rindler coordinates (t, x, y, z) as a function 
of Minkowski coordinates (i, x, y, z) are 

at = e az sinh(ai), az = e az cosh(ai), x = x, y = y (2) 

for region I, and 

at = —e az sinh(at), az — —e az cosh(at), x = x, y = y 

(3) 

for region IV. As we can see from Fig. [TJ there are 
two more regions labeled II and III. To map them we 
would need to switch cosh <->• sinh in Eqns. ([2]) and (|3|). 
In these regions i is a spacclike coordinate and z is a 
timclikc coordinate. However, the solutions of the Klcin- 
Gordon/Dirac equation in such regions are not required 
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to discuss entanglement between the inertial observer Al- 
ice and the accelerated observer Rob. This is so because 
Rob would be constrained to either region I or IV, hav- 
ing no possible access to the opposite regions as they are 
causally disconnected [3, [E l23Tj25| . 

The Rindler coordinates z,t go from — oo to oo inde- 
pendently in regions I and IV. Therefore, each region 
admits a separate quantization procedure with their cor- 
responding positive and negative energy solutions of the 
Klein-Gordon (or Dirac) equations. 



The states 11/ 



,t 



|0) M are free massless scalar 



J-cj/m — u uj,M l u /M 

field modes, in other words, solutions of positive fre- 
quency w (with respect to the Minkowski timclike Killing 
vector df) of the free Klein-Gordon equation: 



|1w)m — u Cj (x 
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|lt2i 2 ) , 



(4) 



where only the time dependence has been made explicit. 
The label M just means that these states are expressed 
in the Minkowskian Fock space basis. 

An accelerated observer can also define his vacuum and 
excited states of the field. Actually, there are two natu- 
ral vacuum states associated with the positive frequency 
modes in regions I and IV of Rindler space-time. These 
are |0)j and |0) IV , and subsequently we can define the 
field excitations using Rindler coordinates (x, t) as 
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(5) 



These modes are related by a space-time reflection and 
only have support in regions I and IV of the Rindler 
space-time respectively. 

However, these Rindler modes are not independent 
of the Minkowskian modes. Indeed we can expand the 
field in terms of Minkowski modes and, independently, in 
terms of Rindler modes. Therefore, the Minkowskian and 
the Rindler set of modes are related by a change of basis 
[TEl [25j . The relationship between the Minkowski Fock 
basis and the two Rindler Fock bases comes through the 
Bogoliubov coefficients, obtained by equating the field 
expansion in Minkowskian modes with the field expan- 
sion in Rindler modes. In general we have that 



"(aK,+^LT* + ««<+^<)- ( 6 ) 

R aR 



The Bogoliubov coefficient matrices a** , (where R = 
I, IV) are given by the Klein-Gordon scalar product be- 
tween both sets of modes 



(7) 



The relationship between modes also establishes a re- 
lationship between the Minkowski annihilation operator 



and the particle operators in Rindler regions I and IV: 



a&j ,m : 
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a aJi,IV + ^ij a wi,l) 



(8) 

On the other hand there exist an infinite number of 
orthonormal bases that define the same vacuum state, 
namely the Minkowski vacuum |0) M , which can be used 
to expand the solutions of the Klein-Gordon equation. 
More explicitly, since the modes n^. have positive fre- 
quency, any complete set made out of independent linear 
combinations of these modes only (without including the 
negative frequency ones u^f*) will define the same vac- 
uum |0) M . 

Specifically, as described in e.g. Refs. |24| - l26| and ex- 
plicitly constructed below, there exists an orthonormal 
basis ^'i^f} determined by certain linear combina- 
tions of monochromatic positive frequency modes, u^f. 
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such that the Bogoliubov coefficients that relate this ba- 
sis {ip^jtp™} and the Rindler basis {u^^u™.} have the 
following form: 



/%• = -(<,<*) =o, 
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(10) 



and analogously for a^ 1 and /3,^' IV interchanging the 
labels I and IV in the formulas above. In this expressions 



tanh 7%, 



exp(— Tru)i/a), 



(11) 



and the label s in r Sj i has been introduced to indicate that 
we are dealing with a scalar field. In this expression and 
in what follows, we will use Planck units (h = c = G = 1). 

In this fashion a mode tp^. (or a mode ip 1 ^ 1 ) expands 
only in terms of mode of frequency u)j in Rindler regions 
I and IV and for this reason we have labeled ip™. and 
ip'^f 1 with the frequency ujj of the corresponding Rindler 
modes. In other words, we can express a given monochro- 
matic Rindler mode of frequency u>j as a linear superpo- 
sition of the single Minkowski modes and ip'^* or as 
a polychromatic combination of the positive frequency 
Minkowski modes vJ^. and their conjugates. 

Let us denote a w and aj, the annihilation and creation 
operators associated with modes ijj^ (analogously we de- 
note a' u . and aj. the ones associated with modes ip'^f)- 
The Minkowski vacuum |0) M , which is annihilated by all 
the Minkowskian operators (Z&^m, is also annihilated by 
all the operators a UJj and a'^ , as we already mentioned. 
This comes out because any combination of Minkowski 
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annihilation operators annihilates the Minkowskian vac- 
uum. 

Due to the Bogoliubov relationships (jTU)) being diago- 
nal, each annihilation operator a Ui can be expressed as 
a combination of Rindlcr particle operators of only one 
Rindlcr frequency w,: 

a Ui = coshr M a Wi I - sinhr s ,i a* . jIV , (12) 

and analogously for a' u . interchanging the labels I and 
IV. 

An analogous procedure can be carried out for 
fermionic fields (e.g. Dirac fields). We can use linear 
combinations of monochromatic solutions of the Dirac 
equation ip™ a and ^™ g (and their primed versions) 
built in the same fashion as for scalar fields: 
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where u^. a and a are respectively monochromatic so- 
lutions of positive (particle) and negative (antiparticlc) 
frequency ±(Dj of the massless Dirac equation with re- 
spect to the Minkowski Killing time. The label a ac- 
counts for the possible spin degree of freedom of the 
fermionic field 1 . 

The coefficients of these combinations are such that for 
the modes ip Ui ,a and rpu ita ^ ne annihilation operators are 
related with the Rindlcr ones by means of the following 
Bogoliubov transformations [H, ]2]i [H| : 

Coj,,a = cosrd,iCi, Wj , CT - sinrd 1 i4v> i ,-a> 

dli,(7 = C0S r d,i 4v.Wi.o- + silir 'd,iCl jWij _ cr , (14) 

and analogously for cf u . and d'£. interchanging the labels 
I and IV, where 

tanrd.i = cxp(— Ttcji/a). (15) 

Here c Uit<T , <i WijCr represent the annihilation operators of 
modes ip^. a and t/>^ a for particles and antiparticles re- 
spectively. The label d in ra.-s has been introduced to 
indicate Dirac field. The specific form for i/j^. a and 
"ip^. „ as a linear combination of monochromatic solu- 
tions of the Dirac equation can be seen, for instance, in 
[13, HH among many other references. Notice again that, 
although we are denotating a Ui , c Ui>a , d UiiCT the operators 
associated with Minkowskian modes, those modes are not 
monochromatic, but a linear combination of monochro- 
matic modes given by © and (|13l) . 



As we are going to discuss fundamental issues and not 
an specific experiment, there is no reason to adhere to 
a specific basis. Specifically, if we work in the bases 
(f^|) . (|T5|) for Minkowskian modes we do not need to carry 
out the single mode approximation 0, [l2| in which one 
single mode of Minkowski frequency Wj was expressed as 
a monochromatic combination of Rindlcr modes of the 
same frequency. This approximation has allowed pioneer- 
ing studies of correlations with non inertial observers, but 
it is based on misleading assumptions on the characteris- 
tics of Rob's detector, being partially flawed. In any case, 
discussing the validity of the single mode approximation 
is not the aim of this work. A complete discussion of the 
problems associated with the single mode approximation 
and how to overcome them is in course of completion 
and will be reported elsewhere [29|. As far as this work 
is concerned it is enough to say that we are not using any 
similar approximation. 

A. Vacuum and first excitation for a scalar field 

The Minkowski vacuum state of the field |0) M is anni- 
hilated by the annihilation operators a,& it M as well as by 
the operators a UJi and a' u .. For the excited states of the 
field, we will work with the orthonormal basis {ip^ , fp'^f} 
defined in ^ such that 

IW M = <|0) M , |1U)m = <|0) m (16) 

are solutions of the free Klein-Gordon equation which are 
not monochromatic, but linear superpositions of plane 
waves of positive frequency ljj . 

As shown in (241 . |25| , we can express the Minkowski 
vacuum state in terms of the Rindlcr Fock space basis, 
|0} M = <8>i l<U) M , where 

l<U) M = coghr , ^(tanhTv,)™ \n^) l |n Ui ) IV . (17) 

s,z n— 

It is straightforward to check that this vacuum is, indeed, 
annihilated by the operators a u>i and a! u)i . 

The Minkowskian one particle state |l Wi ) M (in the ba- 
sis (|9|) ) results from applying the creation operator aj, . to 
the vacuum state. We can also translate it to the Rindlcr 
basis 

1 

IWm " (coshr s , !; ) 2 

OO 

x ^(tanhr s ,i)Vn + 1 \n + |n Wj ) IV . (18) 

n=0 

The mode |1^.)m is analogous but swapping the labels I 
and IV. 



B. Vacuum and first excitation for a Dirac field 

1 Throughout this work we will consider that the spin of each mode 

is in the acceleration direction and, hence, spin will not undergo F °r simplicity in the notation we are going to present 

Thomas precession due to instant Wigner rotations @, [27| . Only the construction of the modes if}^. and tfj^. in @ 
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since, as we will show later when we build the entangled 
state, they are the only ones of importance in this analy- 
sis. In any case, the construction of the primed modes for 
the Dirac field is analogous, but we would need to expand 
the notation below to include the possibility of antipar- 
ticles in Rindler region I and particles in Rindler region 
IV being careful with all the anticommutation subtleties 
typical for fcrmionic fields. 

As for the scalar case, the vacuum state of the field 
|0} M is annihilated by the annihilation operators o-,M 
and (i&^cr.M for all uii,u as well as by the operators c Ui!<7 
and dcj^a- for all u)i,a. 

For the excited states of the field, we will work with 
the orthonormal basis (fT3j) such that 

K>M =C Lr|0> M (19) 

are positive frequency solutions of the free Dirac equation 
which are not monochromatic, but linear superpositions 
of plane waves of positive frequency a)j. 

Let us introduce some notation for the Rindler field 
excitations that will follow the same convention as in 

EES mi. 
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~~ L I,w«,o 
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\(T Wi ) 
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- d f 

~ "iV.Wj 


,J0) IV , 




\Pu 


Jl 


- C f 


cU,JO)l = -cU,|4a>, 


, t l°>i> 


\P<*i) 


IV 


- d f 




,4.^iv,ui,t 0)iv 



(20) 



where p Ui represents the spin pair state in the mode with 
frequency Wj. Notice that due to the operator anticom- 
mutation relations, 

\<Tu t )l K,)lV = t iu i ,<r4v*, i ,*' l°)lV 

= -4v, WjlCT '4, Wi ,J0)il0)iv. 

4v, Ui ,a' K>I |0)l V = - K>I |Oiv • ( 21 ) 

As it can be seen in [15[, the projection onto the un- 
primed sector of the basis (fT3")) of the Minkowski vacuum 
state written in the Rindler basis, is as follows 

|0 Wi ) M = (cosr d) i) 2 |0) I |0) IV 

+ sinr d ,i cosr dil (|t^)i |Wrv + l-Wi IWrv) 
+ (smr d ,*) 2 |p Wi } I |p Wi ) Iv . (22) 

It is straightforward to check that the vacuum is annihi- 
lated by Cui ua and d Uit(7 simply using (fl4|) and applying 
both operators to (j2"2")l . 

The one particle state (projected onto the sector ip^}. 
of (1131) ) in the Rindler basis can be readily obtained by 
applying the particle creation operator cj,. a to (j2"2")l : 

IW M = cosr cU |t«*)i |0) IV + sinr d ,i \p Ui ) 1 |tw 4 )rv ' 
I-Wm = cosr d.» |0) IV - sinr d . 4 (^Jj || Wl ) IV . 

(23) 



C. Entanglement degradation due to Unruh effect 

We will now summarize the results that have been ob- 
tained concerning the effects of an uniform acceleration 
on quantum correlations. 

Let us first consider the following maximally entangled 
state for a scalar field: 

|*> B = i(|0) A |0) R + |l) A |l„ a ) R ), (24) 

where the label A denotes Alice's subsystem and R de- 
notes Rob's subsystem. In this expression, |0) AR repre- 
sents the Minkowski vacuum for Alice and Rob, |1) A is an 
arbitrary one particle state excited from the Minkowski 
vacuum for Alice, and the one particle state for Rob is 
expressed in the basis © and characterized by the fre- 
quency ojr observed by Rob. 

The election of the modes |1 WR ) instead of |1^ R ) to 
build the maximally entangled state is not relevant since 
choosing the primed modes would just be equivalent to 
say that Rob is in region IV instead of in region I, and 
there is complete symmetry in all the analysis for both 
cases. 

Since the second partner (Rob) — who observes the bi- 
partite state (|24p — is accelerated, it is convenient to map 
the second partition of this state into the Rindler Fock 
space basis, which can be computed using equations p7|) 
and (|18[) . and rewrite it in the standard language of rel- 
ativistic quantum information (i.e., naming Alice to the 
Minkowskian observer, Rob to a hypothetic observer in 
Rindlcr's region I and AntiRob to a hypothetic observer 
in Rindler Region IV): 

i*) s = E^^(i°)aKk)hKk)r 

+ ^^|1)aI" + W r Kh)r), (25) 

where 

tanhr s = cxp(— 7rwR,/a). (26) 

The same can be done in the case of a Dirac field. Let 
us now consider the following maximally entangled state 
for a Dirac field in the Minkowskian basis 

|*> d = ^(|0> A |0> R + |t) A |W R )- (27) 

As for the bosonic case, if Rob, who observes this bi- 
partite state, is accelerated, it is convenient to map the 
second partition of this state into the Rindler Fock space 
basis, which can be computed using Eqns. (|2"2")) and (|2"3"|) . 
The explicit form of such state can be seen in [l5j . 

Notice that we have chosen a specific maximally en- 
tangled state (|2"T)) of all the possible choices. This elec- 
tion has no relevance since in (20| it was shown the uni- 
versality of the degradation of fcrmionic entanglement. 
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All fcrmionic maximally entangled states are equally de- 
graded by the Unruh effect, no matter what kind of max- 
imally entangled state is (either occupation number or 
spin Bell state), or even if we work with a Grassmann 
scalar field instead of a Dirac field. 
Let us denote 

ParrH^sX^sI, ^ Rft H*d)(*d|, (28) 

the tripartite density matrices for the bosonic and 
fcrmionic cases, in which we use the Minkowski basis 
for Alice and the Rindler basis for Rob-AntiRob. One 
could ask what is the physical meaning of each of these 
three 'observers'. Alice represents an observer in an in- 
ertial frame. For Alice the states and (|2T|) are maxi- 
mally entangled. Rob represents an accelerated observer 
moving in a x = a -1 trajectory in Region I of Rindler 
space-time (as seen in Fig. [1} who shares a bipartite en- 
tangled state (|24j) or (|27]l with Alice. AntiRob represents 
an observer moving in a x = a" 1 trajectory in Region IV 
with access to the information to which Rob is not able 
to access (at least classically) due to the presence of the 
Rindler horizon. 

In the standard Unruh entanglement degradation sce- 
nario [1, [lH , as Rob is not able to access AntiRob part of 
the system we must trace over AntiRob degrees of free- 
dom when accounting for the quantum state shared by 
Alice and Rob. This provokes, for instance, the observa- 
tion of a thermal bath by Rob while Alice observes the 
Minkowski vacuum as it can be seen elsewhere [H, [TH, [25| . 
As a consequence the state becomes mixed, which causes 
some degree of correlation loss in the system AR as we 
increase the value of the acceleration a. In references 
0, H, [H, [l]l GHim it is studied how this phenomenon 
affects the entanglement for different fields. 

It has been also studied 0, [U [22j the correlation 
trade-off among the all possible bipartitions of the sys- 
tem, namely, Alice- Rob (AR), Alice- AntiRob (AR), and 
Rob-AntiRob (RR). 

Bipartition AR is the most commonly considered in the 
literature. It represents the system formed by an inertial 
observer and the modes of the field which an accelerated 
observer is able to access. The second bipartition (AR) 
represents the subsystem formed by the inertial observer 
Alice and the modes of the field which Rob is not able 
to access due to the presence of an horizon as he acceler- 
ates. Classical communication between the two partners 
is only allowed for the bipartitions AR and AR. We will 
call these bipartitions 'Classical communication allowed' 
(CCA) from now on. These bipartitions are the only 
ones in which quantum information tasks arc possible to 
be performed. 

On the other hand, no quantum information tasks can 
be performed using RR correlations since classical com- 
munication between Rob and AntiRob is not allowed. 
Anyway, studying this bipartition is still necessary to give 
a complete description of the behavior of the correlation 
created between the space-time regions separated by the 
horizon. 



As it is commonplace in quantum information, the par- 
tial quantum states for each bipartition are obtained by 
tracing over the third subsystem 



p AR 




p AR 


= Tr R/ r 




= Tr A p 



In the cases AR and AR, there are physical arguments 
to justify the need for this 'tracing over' beyond mere 
quantum information considerations, namely, Rob will 
never be able to access region IV of the space-time due 
to the presence of the Rindler horizon so that R (Region 
IV) must be traced out. Likewise, AntiRob is not able 
to access region I because of the horizon and hence R 
(Region I) must be traced out. For the subsystem RR 
this tracing over subsystem A corresponds to the stan- 
dard procedure for analyzing correlations between two 
parts of a multipartite system. The properties of the cor- 
relations among these subsystems has been analyzed in 
the literature, showing a completely different behavior of 
quantum correlations for the CCA bipartitions depend- 
ing on whether the system is fcrmionic or bosonic. 

For fcrmionic fields, quantum correlations are con- 
served as Rob accelerates 

Specifically, as entan- 
glement in the bipartition AR is reduced, entanglement 
in the system AR is increased. In the limit of a — > oo 
some entanglement survives in all the bipartitions of the 
system. 

For the scalar field the situation is radically different, 
namely, no entanglement is created in the CCA biparti- 
tions. Moreover, the entanglement in the AR bipartition 
is very quickly lost as Rob accelerates, even if we artifi- 
cially limit the dimension of the Hilbert space [22j ■ 

This different behavior, thought at the beginning to be 
a consequence of the finite dimensionality of the fcrmionic 
Hilbert space, was demonstrated to be ruled only by 
statistics |20l - l22l |. which plays a crucial role in the phe- 
nomenon of Unruh entanglement degradation. The role 
of statistics is so important that, for fermions, the be- 
havior of quantum correlations has been proven to be 
universal [201 ] . Also, the survival of entanglement for 
the fermionic case, is arguably related to statistical cor- 
relations [10, [2f[. All these aspects will be discussed 
in depth later on, when we present the results for the 
Schwarzschild black hole. 



III. THE "BLACK HOLE LIMIT": 
TRANSLATION RINDLER-KRUSKAL 

In this section we will study a com plet ely new setting 
using the tools learned from (j, H, [2f|. We will prove 
in a constructive way that the entanglement degradation 
in the vicinity of an eternal black hole can be studied 
in detail with these well-known tools. By means of the 
construction shown below we will be able to deal with 
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new problems such as computing entanglement loss be- 
tween a free-falling observer and another one placed at 
fixed distance to the event horizon as a function of the 
distance, studying the behavior of quantum correlations 
in the presence of black holes. We will also show that 
the entanglement loss produced by an eternal black hole 
shows universality. 

To begin this section let us work a little bit with the 
Schwarzschild metric 



ds 2 



2 m 



dt 2 



2 m 



dr 2 +r 2 dn 2 , 



(30) 

where m is the black hole mass and dtt 2 is the line el- 
ement in the unit sphere. Due to the symmetry of the 
problem we are going to restrict the analysis to the radial 
coordinate. To shorten notation let us write the radial 
part of metric as 



ds 2 = -/dr + f~ dr 



(31) 



where / = 1 — 2m/ r. 

We can choose to write the metric in terms of the 
proper time to of an observer placed in r = ro as fol- 
lows, 



d.s 2 = -Zdtg + /"Mr 3 
Jo 



(32) 



where /o = 1 — 2m/ tq. The relationship between <o and t 
is given by the norm of the timclikc Killing vector £ = d t 
in r = r , namely t = \fja~t. 

We can now change the spatial coordinate such that 
the new coordinate vanishes at the Schwarzschild radius 
r = Rg = 2m. Let us define z in the following way 



2m 



z 

8m 



/ 



{Kzf 



1 + [nzf 



(33) 



with k = l/(4m) being the surface gravity of the black 
hole. Then the metric (|3"2"j) results 



d,s 2 = - 



1 (kz) 2 

To i + M 



;dt 2 + [l + (KZ) 2 ] dz 2 



(34) 



Near the event horizon (z w 0), we can expand this met- 
ric to lowest order in z and approximate it by 



d^ 



dtl + dz 2 



(35) 



which is a Rindlcr metric with acceleration parameter 
k/VJq. 

On the other hand, Eq. (|35|) represents the metric 
near the event horizon in terms of the proper time of an 
observer placed at r = ro . The next step is giving a phys- 
ical meaning to this Rindler-like acceleration parameter. 
For this, we need to compute the proper acceleration of 
a Schwarzschild observer placed at r = ro, which is, in- 
deed, different from n (as k would be the acceleration of 



an observer arbitrarily close to the horizon as seen from 
a free- falling frame). 

To compute a for this observer as seen by him- 
self (proper acceleration) we must start from the 
Schwarzschild metric. The value of the proper acceler- 
ation for an accelerated observer at arbitrary fixed po- 
sition r is a = y/aJ/aJ 1 where a/ 1 = v v V '„u M is the ob- 
server 4-acceleration at such position, whereas is his 
4- velocity. 

The 4- velocity for a Schwarzschild observer in an arbi- 
trary position r is 



(36) 



where £ = dt is the Schwarzschild timelike Killing vector. 
As £ M = (1,0,0,0) in Schwarzschild coordinates, then 
If I = - V7, and therefore v" = ^/y/J. Thus, we 

can compute the acceleration 4- vector 



(37) 



Taking into account that £ M is a Killing vector and, there- 
fore, it satisfies V M £„ + V„£ M = 0, we easily obtain 



ld^\ 2 d^f 1 



(38) 



Hence, since g rr = /, the proper acceleration for this 
observer is 



W) 5 



4/ 



For an observer placed at r = ro 

K 



a 



(Wo)* 



(39) 



(40) 



We know from Q that 1 - / = [1 + (kzq) 2 }- 1 . So, 
if the observer in r = ro is close to the event horizon 
(ro « Rs), then, to lowest order, 1 + (kzq) 2 rj 1 and 



(41) 



«/v/o 



Therefore, under this approximation, we can re-write 
(1351) as 



ds 2 



(a z) dtl + dz 2 . 



(42) 



This shows that the Schwarzschild metric can be ap- 
proximated, in the proximities of the event horizon, by 
a Rindler metric whose acceleration parameter is the 
proper acceleration of an observer resisting in a position 
ro close enough to the event horizon. 
This approximation holds if 



(A) « 1 



(43) 
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Figure 2: Kruskal space-time conformal diagram showing tra- 
jectories for Alice, Rob and AntiRob. i° denotes the spatial 
infinities, i~ , i + are respectively the timelike past and future 
infinities, and ^ + are the null past and future infinities 
respectively, and M'^ are the event horizons. 



<9 t - oc (d u + d v ). The parameter t for this timclikc 
vector corresponds to the proper time of a free- 
falling observer close to the horizon, and it is anal- 
ogous to the Minkowskian timelike Killing vector. 
Positive frequency modes associated to this time- 
like vector define a vacuum state known as the 
Hartlc-Hawking vacuum |0) H , which is analogous 
to |0) M in the Rindler case. 

dt oc (ud u — vd v ). It is the Schwarzschild timelike 
Killing vector, which (when properly normalized) 
corresponds to an observer whose acceleration at 
the horizon equals the surface gravity k of the black 
hole with respect to a Minkowskian observer, or, in 
other words, with proper acceleration eto sa n/y/fo 
close to the horizon. The vacuum state correspond- 
ing to positive frequencies associated to this time- 
like Killing vector is called the Boulwarc vacuum 
|0) B . This state is analogous to the Rindler vac- 
uum |0) T . 



or, in other words, if 



(44) 



where Ao = ro — Rs is the distance from ro to the event 
horizon. In the limit r — > i?s wc obtain that /q — > and, 
from (|41[) . ao — > oo. This shows rigorously that being 
very close to the event horizon of a Schwarzschild black 
hole can be very well approximated by the infinite accel- 
eration Rindler case, as it was suggested in pi HL l20l l2l| . 
This also enables us to study what would happen with 
the entanglement between observers placed at different 
distances of the event horizon as far as the Rindler ap- 
proximation holds. 

Now let us identify again who is who in this new sce- 
nario. For this, we introduce the null Kruskal-Szeckeres 
coordinates 

u = -k' 1 exp[-n(t-r*)}, v = k' 1 exp[K,(t+r*)}, (45) 

where r* = r + 2mlog|l — r/2m\. In terms of these 
coordinates the radial part of the Schwarzschild metric is 



ds 2 = — e- 2Kr dudv, 
2nr 



(46) 



where r is implicitly defined by (|45[) . The Penrose di- 
agram for this maximal analytic extension is shown in 
Fig. [2] In this coordinates, near the horizon the metric 
can be written to lowest order as 



ds 2 



-e 1 dudv 



(47) 



and uv = —(kz) . 

Hence, there are three regions in which we can clearly 
define physical timelike vectors respect to which we can 
classify positive and negative frequencies: 



• There is another timelike Killing vector —dt (as in 
Rindler) for region IV that will allow us to define 
another Boulware vacuum in region IV. We will call 
it AntiBoulwarc vacuum |0)g, analogous to |0) IV in 



the Rindler case. 



Now, in this scenario, |l<a) H = H |0) H are free scalar 
field modes, in other words, solutions of positive fre- 
quency Co with respect to 9 t - of the free Klein-Gordon 
equation close to the horizon 



I 1 ^)h = u * a 



f 



>2Cj 



(48) 



The label H just means that those states are expressed 
in the Hartle-Hawking Fock space basis. 

An observer located at a fixed distance from the black 
hole can also define his own vacuum and excited states 
of frequency u> respect to the Killing vector dt- Actu- 
ally, there are two natural vacuum states associated with 
the positive frequency modes in both sides of the hori- 
zon these are |0) B and |0) B , vacua for the positive fre- 
quency modes in regions I and IV respectively (Fig. [2|. 
Subsequently, for a scalar field, we can define the field 
excitations as 



f 



l 1 ">B = a L,B l°>B = U u> « 



I 1 ")b = a l,g I°)b = U u « ~^f' 



(49) 



Then, the analogy between the Rindlcr-Minkowski and 
the Boulwarc-Hartlc-Hawking states, and their relation 
with the standard Alice-Rob- AntiRob notation is as fol- 
lows: 



|o) R 
|o)r 

|0)a 



|0) IV |0)g 



/B ' 



B • 



|0) M |o) H - 



(50) 
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The change of basis between Hartlc-Hawking modes and 
Boulware modes is completely analogous to the change 
of basis between Minkowskian modes and Rindler modes 
with an acceleration parameter do = k/ vTo- 

In the same fashion as for Rindler we define an or- 
thonormal basis of Hartlc-Hawking scalar field modes 
{V'Sj ) V^} whose elements are superpositions of positive- 
frequency solutions Uq. of the Klein-Gordon equation 
with respect to the Kruskal time t such that each element 
corresponds to Boulware modes of one single frequency 
in the Kruskal regions I and IV (u®. and it B *). The same 
can be done for the Dirac field. 

We can express the Hartlc-Hawking vacuum state in 
terms of the Boulware Fock space basis. To do so we 
use what we learned from the Rindler case. Taking into 
account that |0) H = 0^ ICL^jj, we have that 

1 °° 

I°^)h = r ^(tanhg s . 4 )" \n Ui ) B \n Ui ) B , (51) 

cosng Sjl n=Q 

where 

tanhq M = exp {^-Tt^ffo Wj/ftJ ■ (52) 

The unprimed Hartle-Hawking one particle state in the 
basis {ip^ results from applying the correspond- 
ing creation operator to the vacuum state. We can also 
translate this state to the Boulware basis: 

11 \ - 1 

(cosh^,:) 2 

CO 

x y^(tanhff s ,i) T Vn+ l\n+ W B |n Wj ) B . (53) 

n=0 



scenario which looks like follows, for fcrmions and bosons, 
in the basis of a free-falling observer (Alice) 

|*) 8 «(|0) a |0) r +|1)J1 Wr )r), (57) 
|f) d =^=(|0) A |0) R +|t) A |W R )- (58) 

This bipartite system consists in two subsystems, the first 
one is going to be observed by Alice, who is free-falling 
into the black hole and close to the event horizon, and 
the second one will be observed by Rob, who is near the 
event horizon at r = r$ w Rg. Therefore, the second 
partner who observes the bipartite states (|57|) and (f5"5)) 
describes them using the Boulware basis, so that it is 
convenient to map the second partition of these states 
into the Boulware Fock space basis. 

Following the notation (|5T)|) , to analyze the correlations 
among the bipartite subsystems we need to trace out the 
third subsystem analogously to what we did in (f29|) : 





Tr a p ARR 


p AR = 


Tr R p ARR 




Tr A p ARR 



It can be seen in Fig. [5]that all the information beyond 
the event horizon cannot be accessed by Rob. Actually, 
what happens beyond the horizon is determined by the 
information that Rob can access along with the informa- 
tion that AntiRob can access. In this context it makes 
sense to say that studying the system p RR gives an idea 
of the correlations across the horizon. 



IV. CORRELATIONS BEHAVIOR 



The Hartle-Hawking vacuum (projected onto the un- 
primed sector) for the Dirac case is expressed in the Boul- 
ware basis as follows 

l<U> H = (cos 9cM ) 2 \0 U ,) B \0 Ui ) B 

+ sin qd,i cos gd,i (ITo^b I4-Wj)g + I-Wb IWb) 
+ (sing di? ) 2 \p Ui ) B \p Ui )a , (54) 

whereas the projected Hartlc-Hawking one particle state 
is expressed in the Boulware basis as 

IW H = COS q d ,i |tw i ) B l°"i)g + sin 9d. l \p Ui ) B |W B > 
|4-wi) H = cos 1d,i \lu>i) B |0 W4 )g - sing dji \p Ui ) B \lu,,) B , 

(55) 

where this time 

tang<j,i = cxp {^-tt^JJu Ui/n) . (56) 

Thus, in this new scenario, we can consider a bipartite 
state analogous to the states (|24p and (f2"T)) for the Rindler 



In this section we will use the machinery we already 
have from the Rindler set-ups to compute the entangle- 
ment degradation as a function of the position of Rob. 

First we will consider that Rob's frequency ur is mea- 
sured in natural units adapted to each black hole. This 
will show how modes of different frequencies suffer dif- 
ferent correlation degradation. It will also show how less 
massive black holes produce a higher degradation than 
the heavier ones. Furthermore, this analysis will show 
the universality of the phenomenon of the Hawking en- 
tanglement degradation for Schwarzschild black holes. 

After that, we will analyze the different degree of en- 
tanglement degradation experimented by an observer of 
fixed Boulware frequency standing at fixed distances 
from the event horizon for different black hole masses. 

In the following subsections we will see that all the 
interesting behavior happens in regions in which the 
Rindler approximation (|4"2"j) is valid. Specifically, in the 
plots below, the values of the distance to the horizon from 
which the interesting entanglement behavior appears are 
Ao < 0.05i?s hi ah the cases considered in this section 
for which, consequently, the approximation (I42[) holds. 
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A. Adapted frequency 

In terms of the mode frequency measured by Rob (writ- 
ten in units natural to the black hole, i.e. in terms 
of the surface gravity k) and his position measured in 
Schwarzschild radii, 



O = 27rcjp{,/ft = 87rmwR, 
Ro = r /Rs = r /(2m), 



Eqns. ([52jl and (|56|) can be written as 



1 ~ R~a 



fi / 1 
tang d = cxp ( - — J 1 - — 



(60) 
(61) 



(62) 
(63) 



showing that the phenomenon of Hawking entangle- 
ment degradation presents universality, which is to say, 
if the frequency is measured in natural units, every 
Schwarzschild black hole behaves in the same way, as 
expected. 



1. Quantum correlations 

We will use the negativity (AT) to account for the quan- 
tum correlations between the different bipartitions of the 
system. It is an entanglement monotone sensitive to dis- 
tillablc entanglement. The negativity is defined as the 
sum of the negative eigenvalues of the partial transpose 
density matrix for the system, which is defined as the 
transpose of only one of the subsystem q-dits in the bi- 
partite density matrix. For a general density matrix of a 
bipartite system AB, 



PAB 



^2 Pv kl 

ijkl 



K)a \j)b Wa (l\ 



B ■ 



the partial transpose is defined as 

PAB = Pv kl I*) A I0 B ( fc U 01b ■ 
ijkl 

If Aj are the eigenvalues of p^B' then 

i Ai<0 



(64) 



(65) 



(66) 



Hence, to compute it, we will need the partial trans- 
pose of the bipartite density matrices (|59l) . The details 
associated to the diagonalization of the partial trans- 
posed density matrices for each subsystem are technically 
very similar to the Rindler case, and are not of much in- 
terest for the purposes of this article. All the technical 



aspects of such calculations can be found in [21( for Dirac 
and scalar fields. The results of those calculations are 
shown in Figs. |3]to|ni In Figs. [3] and 2] we can see the 



behavior of the negativity on the CCA bipartitions for 
different values of Rob's frequency fi. 

For the scalar field we can see that as Rob is closer to 
the event horizon the entanglement shared between Alice 
and Rob decreases. In the limit in which Rob is very close 
to the horizon, entanglement is completely lost. With 
the study performed in this work we can see the func- 
tional dependence of the entanglement with the distance 
to the horizon. As seen in the figures, the degradation 
phenomenon occurs in a narrow region very close to the 
event horizon. If Rob is far enough from the black hole he 
will not appreciate any entanglement degradation effects 
unless either the mass of the black hole or the frequency 
of the mode considered are extremely small. There must 
be, indeed, a minimum residual effect associated to the 
Hawking thermal bath experienced in the asymptotically 
flat region of the spacetime, far from the region in which 
this approximation is valid, but it is unnoticeable small. 
Certainly, as it will be seen in fig. [9] and the discussion 
below, even very close to the horizon no effective en- 
tanglement degradation occurs for physically meaningful 
values of mass and frequency. 

If we keep the frequency measured by Rob cjr con- 
stant, n will grow proportional to the black hole mass. 
With this in mind, Fig. [3] shows that the degradation 
is stronger for less massive black holes. This result is 
consistent with the fact that the Hawking temperature 
increases as the mass of the black hole goes to zero. In 
the next section (specifically in Fig. [9]) wc will show that 
this is not an effect of choosing natural units, when an 
observer is at a fixed distance of a black hole, the degra- 
dation will be higher for less massive black holes. 

In any case, for the scalar field, the entanglement in 
the system AR is completely degraded when one of the 
observers is resisting very close to the event horizon of 
the black hole. Hence, in this scenario, no quantum in- 
formation resources can be used (for instance to perform 
quantum teleportation or quantum computing) between 
a free-falling observer and an observer arbitrarily close 
to an event horizon. Moreover, no entanglement of any 
kind is created among the CCA bipartitions of the system 
(the ones who can classically communicate). Therefore, 
all useful quantum correlations between a free-falling ob- 
server and an observer at the event horizon are lost due 
to the Hawking effect degrading all the entanglement in 
the system. 

For the Dirac field (Fig. 3} something very different 
happens. We see that correlations in the bipartition AR 
decrease to a certain finite limit, which means that there 
is entanglement survival even when Rob is asymptoti- 
cally close to the event horizon. This survival is a well 
known phenomenon in the Rindler case 

dm- At th ° 

same time that entanglement is destroyed in the AR bi- 
partition, entanglement is created in the complementary 
AR bipartition so that negativity in the CCA biparti- 
tions fulfills a conservation law regardless of the distance 
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Position of Rob (ro/i? s ) 




1.05 



1.04 1.03 1.02 1.01 

Position of Rob (rg/i? s ) 



Figure 3: Scalar field: Entanglement of the system Alice- Rob 
as a function of the position of Rob for different values of Q. 
Entanglement vanishes as Rob approaches the Schwarzschild 
radius while no entanglement is created between Alice and 
AntiRob. The smaller the value of f2 the more degradation is 
produced by the black hole. 



Figure 4: Dirac field: Entanglement Alice-Rob (blue solid 
line) and Alice- AntiRob (red dashed line). Universal conser- 
vation law for fermions is shown for different values of SI. The 
entanglement degradation in AR is quicker when Q, is smaller. 
The maximum degradation is not maximal and its value is in- 
dependent of fi. 



to the event horizon and the mass of the black hole 

A^AR+AA A R, = i (67) 

The nature of this entanglement and the survival of cor- 
relations, even in the limit of positions arbitrarily close 
to the horizon, is discussed in [U [22[ for the Rindlcr 
case. When we deal with fcrmionic fields there are cor- 
relations that come from the statistical fermionic nature 
of the field which we cannot get rid of. The hypothe- 
sis is that this entanglement, which is purely statistical, 
is the second quantized version of the statistical entan- 
glement disclosed in [3(j. Here we see that the same 
conclusions drawn in that case can be perfectly applied 
to the Schwarzschild black hole case. 

About the dependence of the entanglement degrada- 
tion on the frequency of the Boulware mode, Fig. [3] 
shows that, for a scalar field, the loss of entanglement be- 
tween a free falling observer and an observer outside but 
very close to the event horizon (AR) is greater for modes 
of lower frequency. This makes sense because, energeti- 
cally speaking, it is cheaper to excite those modes and, 
therefore, they are more sensitive to the Hawking thermal 
noise. For a Dirac field (Fig. @| we see a similar behav- 
ior, namely, lower frequencies are less protected against 
entanglement degradation due to Hawking effect. How- 
ever, the surviving entanglement in the limit in which 
Rob is infinitely close to the event horizon is not sensi- 
tive to the frequency of the mode considered; remarkably, 
the entanglement decays up to the same finite value for 
all modes. This is in line with the idea that the entan- 
glement that survives the event horizon is merely due 



to statistical correlations, and the only information that 
survives when Rob is exactly at the horizon is the fact 
that the field is fermionic as it is discussed in (20l - [22| . 

From Figs. [3] and @] we can also conclude that all 
the relevant entanglement degradation phenomena is pro- 
duced in the proximities of the event horizon so that the 
Rindlcr approximation that we are carrying out is valid 
(Eq. 23| . We can also see that the degradation is small 
even in regions in which the approximation still holds. 
Therefore for longer distances from the horizon the pres- 
ence of event horizons is not expected to perturb entan- 
gled systems. 

In Figs. [5] and |3] we can see the behavior of the nega- 
tivity on the RR bipartition for scalar and Dirac fields re- 
spectively. Here we see that quantum correlations across 
the horizon are created as Rob is standing closer to the 
event horizon. In other words, as Rob is getting closer to 
the event horizon the partial system RR gains quantum 
correlations. This result shows that, when Rob is near 
the horizon, the field states in both sides of the event 
horizon are not completely independent. Instead, they 
get more and more correlated. However, this RR entan- 
glement is useless for quantum information tasks because 
classical communication between both sides of an event 
horizon is forbidden, so no quantum teleportation can be 
performed across the event horizon. It is well known for 
the Rindler case that quantum correlations are created 
between Rob and AntiRob [2l| when the acceleration in- 
creases. Here we see the direct translation to the Kruskal 
scenario. The growth of those correlations encodes infor- 
mation about the dimension of the Fock space for each 
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Figure 5: Scalar field: Entanglement of the system Rob- 
AntiRob (entanglement across the horizon) as a function of 
the position of Rob for different values of Q. Entanglement 
diverges as Rob approaches the Schwarzschild radius. 



Figure 7: Scalar field: Mutual information Alice-Rob (blue 
solid line) and Alice- AntiRob (red dashed line). A conserva- 
tion law (derived from the behavior of purely classical corre- 
lations) is shown. Mutual information AR decreases as Rob 
is closer to the horizon and mutual information AR grows. 



0.6 n = i 




1.05 1.04 1.03 1.02 1.01 1 



Position of Rob (ro/i? s ) 

Figure 6: Dirac field: Entanglement of the system Rob- 
AntiRob (entanglement across the horizon) as a function of 
the position of Rob for different values of Q. Entanglement 
tends to a finite value as Rob approaches the Schwarzschild 
radius. 



field mode [H. 



2. Mutual information 

To account for all the correlations among the different 
bipartitions of the system we will use the mutual infor- 



mation, which accounts for correlations (both quantum 
and classical) between two different parts of a system. 
For a bipartite system AB, it is defined as 

Iab = S A + S B - Sab, (68) 

where Sa, Sb and Sab are respectively the Von Neumann 
entropies S = Tr (p log 2 p) for the individual subsystems 
A and B and for the joint system AB. To compute the 
mutual information for each bipartition we will need the 
eigenvalues of the corresponding density matrices. Again 
the technicalities of this analysis can be found elsewhere 
0, [SI- The results for the CCA bipartitions are shown 
in Figs. [7] and E 

We see here that we obtain the black hole version of the 
mutual information universal conservation law found in 
previous works for the Rindler case plj . Namely, for any 
distance to the horizon or black hole mass it is fulfilled 
that 

^AR + / A r = 2. (69) 

Although, as we can see comparing Figs. [7] and [5J the 
behavior of the mutual information is very similar for 
both fermions and bosons, the origin of this conservation 
law near the event horizon is completely different. 

For scalar fields this conservation near the horizon re- 
sponds to a conservation of classical correlations only. 
This can be deduced from Fig. [JJwhich shows that quan- 
tum correlations drop very quickly as the distance of Rob 
to the horizon decreases and, consequently, the only cor- 
relations left must be classical. However, the conserva- 
tion of classical correlations in the CCA bipartitions has 
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Figure 8: Dirac field: Mutual information Alice-Rob (blue 
solid line) and Alice- AntiRob (red dashed line). A conserva- 
tion law (derived from the behavior of purely quantum corre- 
lations) is shown. Mutual information AR decreases as Rob 
is closer to the horizon and mutual information AR grows. 



Figure 9: Scalar field: Entanglement Alice-Rob when Rob 
stands at a distance of 1 cm and 10 cm from the event hori- 
zon for a fixed frequency o;r = 1.5 Mhz as a function of the 
black hole mass. Notice that, for these values of Aq, the 
approximation holds perfectly for any mass m > 10 solar 
masses. 



to do with the infiniteness of the dimension of the Hilbcrt 
space, as it is shown in p2j . If the dimension of a bosonic 
field is limited to a finite value, classical correlations also 
drop as Rob is closer to the horizon (as quantum corre- 
lations do). 

On the other hand, a Dirac field has a built-in di- 
mensional limit for the Hilbert space of each mode im- 
posed by Pauli exclusion principle. Although previous 
works demonstrated that this limit in the dimension has 
not hing; to do with the behavior of quantum correlations 
IlSmfj it does limit the creation of classical correlations. 
Analogously to what is discussed in [22|, the origin for 
the conservation law (|69p in the fermionic case is a di- 
rect consequence of the quantum correlations conserva- 
tion law ((57]). 

The conclusion is that, although (|6T)]) is universal for 
scalar and Dirac fields in the proximity of an eternal black 
hole, its origin is completely different. For scalar fields it 
responds to a conservation of classical correlations while 
for Dirac fields it is reflecting the quantum correlations 
conservation (pjT)) . 

Mutual information for the RR bipartition does not 
add any new result as it inherits the quantum correlations 
behavior showed in Figs. [3]and|4j 



B. Entanglement degradation dependence on the 
black hole mass 

In this section we will analyze the entanglement degra- 
dation for an observer with the same characteristics in the 
presence of different black holes. To do so we are going 



to use the full dimensional quantities ur and Ao- 

We will consider that Rob's mode frequency is = 
1.5 Mhz, and he is standing at a distance Ao = 1 cm and 
Ao = 10 cm from the event horizon of black holes with 
different masses, while he shares an entangled state ([ST)) 
or (|58[) with a free-falling observer Alice. 

The quantum correlations that Rob and Alice share 
are shown in Figs. [5] and [TU] for scalar and Dirac fields, 
respectively. From these figures we see that for a really 
close distance from the event horizon, only small black 
holes would produce significant entanglement degrada- 
tion. Actually, the degradation decreases very quickly as 
the black hole mass is increased. 

Furthermore, we can see that the effects on the en- 
tanglement decrease very quickly as the distance to the- 
event horizon is increased. This shows that quantum in- 
formation tasks can be safely performed in universes that 
present event horizons since only in the closest vicinity of 
the less massive black holes the Hawking effect impedes 
the application of quantum information protocols. 



V. LOCALIZATION OF THE STATES 

Along this work we have used a plane-wave-like basis 
to express the quantum state of the field for the inertial 
an accelerated observers. These plane wave modes are 
completely delocalized, and therefore, they are not the 
most natural election of modes if we want to think of 
the observers Alice and Rob spatially localized to some 
degree. 
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Figure 10: Dirac field: Entanglement AR (blue continuous 
line) and AR (red dashed line) when Rob stands at a distance 
of 1 cm and 10 cm from the event horizon for a fixed frequency 
ujr = 1.5 Mhz as a function of the black hole mass. Notice 
that, for these values of Ao, the approximation holds perfectly 
for any mass m > 10~ 5 solar masses. 



However, a very similar analysis to the one carried out 
in sections [TT| and IIIII can be performed using a complete 
set of wave packet modes for both the Minkowski and 
Rindlcr solutions of the wave equation. These modes 
can be spatially localized and provide a clearer physical 
interpretation for Alice and Rob, which will eventually 
have to carry out measurements on the field. The way to 
build these wave packet modes can be found elsewhere 

mmm. 

The elements of this basis are defined as a function of 
the plane wave modes (|4]) as 



u Cj,1 



1 



dv e 



-ivl M 



(70) 



where 0) and I label each wave packet. 

We can define creation an annihilation operators asso- 
ciated to these wavepackets aci,i,u> a>t, i m such that 



annihilates the Minkowski vacuum and al , M |0) M = 
|l[i,i) M represents a wavepacket peaked for a frequency 
w and whose spatial localization can be associated to the 
maximum of u^, as a function of x and t. 

A similar analysis can be done for the Rindlcr basis 



.IV 



1 



uj-\-e 



uj-\-e 



dve~ ivl u]7, 



(71) 



u> and I' label each wave packet. We can define creation 
an annihilation operators associated to these wavepack- 
ets a u ,i',R, r (where R = I, IV) such that a u ,z',ij 



annihilates the region Rindler region R vacuum and 

a L v r \ ®)r = \^-u,1')r represents a wavepacket peaked for 
a Rindlcr frequency u> and whose spatial localization can 
be associated to the maximum of v as a function of x 
and t. 

We can compute then the Bogoliuvob transformation 
between the Minkowski wavepackets and the Rindlcr 
wavepackets [32j 



a.u>,i,M 



z ui,l' ,Ci,l a u,l',I 



■p, 



u,l',Ci,l %,!',IV 



(72) 



Where the Bogoliuvob coefficients are computed in the 
same fashion as for the plane wave case 



a. 



I 

Ul,l' ,ui 7 l 



,.M 



uj,V ,Cj,1 



[*\ 

,l) 



(73) 

It is shown in [32j that, apart from an irrelevant phase 
factor, the Bogoliuvob coefficients are related with (fT0|) 
as follows 



(74) 



that G a (ui,l,u>,l') w S UUa 5u a and 
s^zifl) where l a = l a (l') and uj a 



1 1 is shown in [3j 
uj a (u),l'). 

The key feature of this transformations is that they 
have again a diagonal form. As it can be read from (|74p , 
a Minkowski wavepacket |l W; z') M is connected with a pair 
of Rindlcr wavepackets in regions I and IV. Moreover, the 
functional form of the dependence of this coefficients with 
the acceleration is effectively the same. This analysis 
made for the Rindler and Minkowskian modes can be 
straightforwardly translated to the Boulware and Hartle- 
Hawking modes. A completely analogous analysis can be 
done for the fermionic case. 

Consequently all the conclusions extracted in this ar- 
ticle for delocalised modes are also valid for the localized 
modes defined above. 



VI. CONCLUSIONS 

We have analyzed the entanglement degradation pro- 
duced in the vicinity of a Schwarzschild black hole. 

With this aim, we have carried out a detailed study of 
the Schwarzschild metric in the proximity of the horizon, 
showing how we can adapt the tools developed in the 
study of the entanglement degradation for uniformly ac- 
celerated observers 0, H, [H, [H| to the black hole case. In 
particular, we have shown that, regarding entanglement 
degradation effects, the Rindler limit of infinite acceler- 
ation reproduces a black hole scenario in which Rob is 
arbitrarily close to the event horizon. More importantly, 
we have shown the fine structure of this limit, making 
explicit the dependence of the entanglement degradation 
phenomena on the distance to the horizon, the mass of 
the black hole, and the Boulware frequency ur of the en- 
tangled mode under consideration, while keeping control 
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of the approximation to make sure that the toolbox de- 
veloped for the Rindler case can be still rigorously used 
here. 

By means of this analysis we have seen that all the 
interesting entanglement degradation phenomena due to 
the Hawking effect are produced very close to the event 
horizon of the Schwarzschild black hole. The entangle- 
ment degradation introduced by the Hawking effect be- 
comes quickly negligible as Rob is further away from 
the event horizon. In other words, quantum information 
tasks done far away from event horizons are not per- 
turbed by the existence of such horizons. 

We have also shown that for a fixed Rob's mode fre- 
quency and at a fixed distance from the event horizon 
the entanglement degradation is greater for less massive 
black holes. This is consistent with the fact that the 
Hawking temperature is higher for less massive black 
holes. Furthermore, the Hawking entanglement degra- 
dation is a universal phenomenon in the sense that the 
degradation depends only on Rob's frequency and his 
distance to the horizon in units natural to the black 
hole (namely, the surface gravity for frequencies and the 
Schwarzschild radius for distances). In these units, there 
is no extra dependence on the black hole mass, as ex- 
pected. 

We have been able to adapt all the conclusions drawn 
for the Rindler case to the Schwarzschild scenario. In 
particular, we have seen that bosonic and fermionic en- 
tanglement behave in a very different way in the prox- 
imity of a black hole. As it was known for the Rindler 
entanglement on the CCA bipartitions is com- 
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case 

pletely lost for the scalar field while there is a quantum 
correlation conservation law for the Dirac field. 

In [2(| it was shown that for two different kinds of 
fermionic fields (Dirac fields or Grassmann scalars) and 
also for different maximally entangled states (occupation 
number or spin Bell states) the entanglement in the CCA 
bipartitions behaves exactly the same way. This fact was 
used to argue that it is statistics and not dimensional- 
ity that determines the behavior of correlations in the 
CCA bipartitions in the case of uniformly accelerated 
observers. This study proves that this argument is also 
valid for Schwarzschild black holes, not only in the limit 
in which Rob is on the event horizon, but in the whole 
region in which the interesting entanglement degrada- 
tion phenomena are produced. Therefore, the universal 
fermionic entanglement behavior is also manifest in the 
presence of a black hole. 

For the Schwarzschild case, there also appears the uni- 
versal mutual information conservation law found for 
both scalar and Dirac fields in the Rindler case [2l|. In 
the fermionic case, it is due to a conservation of quantum 



correlations while, for bosons, it only reflects the conser- 
vation of classical correlations that happens in the case 
of infinite dimensional Hilbert spaces for each mode. 

Moreover, as Rob is getting closer to the event hori- 
zon, quantum correlations between modes on both sides 
of the event horizon are created, namely the correlations 
between field modes in region I and IV of the Kruskal 
space-time grow up to a value determined by the dimen- 
sion of the Hilbert space of each mode, which is finite for 
the fermionic case and infinite for the scalar field. 

As discussed previously for the Rindler scenario 
[22|| , and unveiled here for the Schwarzschild black hole 
case, under the hypothesis that the fermionic entangle- 
ment which survives the event horizon is of statistic na- 
ture as in [30j ] . it would not contain any useful informa- 
tion. Therefore if an entangled pair is created close to 
the event horizon (for instance particle/antiparticle cre- 
ation) and one of the subsystems falls into the black hole 
while the other resists close to the horizon, no other en- 
tanglement except for the mere statistical would survive 
the degradation provoked by the Hawking effect. 

The problem of the localization of the Rindler and 
Minkowski modes has also been analyzed, showing that 
the results obtained here can be extrapolated to the case 
in which we consider a complete set of localized wave 
packets as a basis of the Fock space for the inertial and 
accelerated observers. 

The scenario that we have discussed here is that of a 
static eternal black hole in which no dynamics is present. 
The analysis of entanglement degradation due to the dy- 
namical creation of an event horizon in a gravitational 
collapse scenario is under current development and will 
be reported elsewhere. 
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